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Background:  To  improve  our  ability 
to  identify  physiologic  deterioration  caused 
by  critical  injury,  we  applied  nonlinear 
analysis  to  the  R-to-R  interval  (RRI)  of  the 
electrocardiogram  of  prehospital  trauma 
patients. 

Methods:  Ectopy  -free,  800-beat  sec¬ 
tions  of  electrocardiogram  from  31  patients 
were  identified.  Twenty  patients  survived 
(S)  and  11  died  (NonS)  after  hospital  admis¬ 
sion.  Demographic  data,  heart  rate,  blood 
pressure,  field  Glasgow  Coma  Scale  (GCS) 
score,  and  survival  times  were  recorded. 
RRI  complexity  was  assessed  via  nonlinear 
statistics,  which  quantify  entropy  or  fractal 
properties. 


BeSUltS:  Age  and  Held  heart  rate  and 
blood  pressure  were  not  different  between 
groups.  Mean  survival  time  (NonS)  was  129 
hours  ±  62  hours.  NonS  had  a  lower  GCS 
score  (8.6  ±  1.7  vs.  13.2  ±  0.8,  p  <  0.05). 
RRI  approximate  entropy  (ApEn;  0.87  ± 
0.06  vs.  1.09  ±  0.07,  p  <  0.01),  sample  en¬ 
tropy  (SampEn;  0.80  ±  0.08  vs.  1.10  ±  0.05, 
p  <  0.01)  and  fractal  dimension  by  disper¬ 
sion  analysis  (1.08  ±  0.02  vs.  1.13  ±  0.01, 
p  <  0.05)  were  lower  in  NonS.  Distribution 
of  symbol  2  (Dis_2),  a  symbol-dynamics 
measure  of  RRI  distribution,  was  higher  in 
NonS  (292.6  ±  34.4  vs.  222  ±  21.3,  p  < 
0.10).  For  RRI  data,  logistic  regression  anal¬ 
ysis  revealed  ApEn  and  Dis_2  as  indepen¬ 


dent  predictors  of  mortality  (area  under  the 
receiver-operating  characteristic  curve  = 
0.96).  When  GCSMOXOR  was  considered,  it 
replaced  Dis_2  whereas  ApEn  was  retained 
(area  under  curve  =  0.92).  When  Injury 
Severity  Score  was  considered,  it  replaced 
GCSMotor;  ApEn  was  retained. 

Conclusions:  Prehospital  loss  of  RRI 
complexity,  as  evidenced  by  decreased  en¬ 
tropy,  was  associated  with  mortality  in 
trauma  patients  independent  of  GCS 
score  or  Injury  Severity  Score. 
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Remote  triage,  or  the  early  categorization  from  a  distance 
of  patients  according  to  injury  severity,  may  enable 
decision  support  and  increase  situational  awareness  on 
the  future  battlefield.  Conventional  vital  signs  such  as  heart 
rate  and  blood  pressure,  however,  are  not  ideal  for  this  pur¬ 
pose  as  they  have  limited  predictive  accuracy.  Thus,  improve¬ 
ment  in  prehospital  triage  will  require  new  analytical 
approaches.1 
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To  generate  new,  more  informative  vital  signs,  analysis 
of  the  electrocardiogram  (ECG)  has  been  pursued  using  a 
variety  of  techniques.  One  of  these  techniques  (spectral  anal¬ 
ysis)  measures  the  strength  of  periodic  oscillations  in  the 
heart  rate.  Heart-rate  variability  (HRV),  that  is,  variability  in 
the  length  of  the  R-to-R  interval  (RRI)  of  the  ECG,  is  affected 
by  autonomic  nervous  system  activity.2-4  RRI  oscillations  in 
the  high-frequency  range  (those  occurring  at  the  respiratory 
rate,  known  as  respiratory  sinus  arrhythmia)  reflect  vagal 
heart-rate  control.  RRI  oscillations  at  the  low-frequency 
range  (about  once  every  10  seconds)  reflect  both  vagal  and 
sympathetic  control  of  the  heart  rate.  By  taking  the  ratio  of 
the  low-frequency-to-high-frequency  powers  of  these  oscil¬ 
lations,  a  measure  of  the  balance  between  sympathetic  and 
vagal  cardiac  activities  can  be  obtained  (“sympathovagal 
balance”).3  As  one  would  expect,  hypovolemia  has  been 
shown  to  cause  a  compensatory  increase  in  sympathovagal 
balance.5-7  On  the  other  hand,  sympathetic  failure,  mani¬ 
fested  by  a  decrease  in  sympathovagal  balance,  was  pre¬ 
dictive  of  death  in  intensive  care  unit  patients8'9  and  in 
prehospital  trauma  patients.10’11 

Rather  than  measure  the  strength  of  the  periodic  oscilla¬ 
tions  in  the  heart  rate,  an  alternative  approach  to  analysis  of 
RRI  variability  measures  the  amount  of  irregularity,  or  com¬ 
plexity,  in  the  heart  rate.  This  irregularity  is  the  result  of  the 
intricate  and  multilayered  control  mechanisms  that  fine-tune 
the  heart  rate.12  Complexity  is  a  feature  of  normal  cardiovas¬ 
cular  regulation,  whereas  loss  of  complexity  or  “decomplexi- 
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RRI  Complexity  and  Mortality 


fication”  of  the  RRI  has  been  shown  in  several  studies  to  be 
a  feature  of  disease13  and  of  impaired  adaptation  to  physio¬ 
logic  stress.14 

In  this  article,  we  use  the  term  “complexity  analysis”  to 
refer  to  a  family  of  tools  more  properly  termed  nonlinear 
statistics.  Nonlinear  behavior  is  the  rule  in  biology.  Biologic 
processes  cannot  be  described  by  analysis  of  the  simple  sums 
of  component  variables,  but  involve  products  and  powers  of 
these  constantly  interacting  variables.  In  the  cardiovascular 
system,  one  consequence  of  this  nonlinearity  is  the  complex 
irregularity  of  the  RRI.  Nonlinear  statistics  are  well  suited  for 
analysis  of  such  signals. 

We  previously  found  in  two  animal  models  of  hemor¬ 
rhagic  shock  that  RRI  complexity,  as  measured  by  the  Ap¬ 
proximate  Entropy  technique,  decreases  during  shock  and  is 
restored  by  fluid  resuscitation.15'16  Similarly,  others  docu¬ 
mented  a  decrease  in  the  complexity  of  the  RRI  in  human 
volunteers  subjected  to  central  hypovolemia  by  means  of 
lower  body  negative  pressure.17  The  objective  of  the  present 
study  was  to  evaluate  the  ability  of  several  such  methods  to 
discriminate  between  survivors  and  nonsurvivors  of  trauma  in 
the  prehospital  setting.  We  hypothesized  that  loss  of  RRI 
complexity  is  associated  with  mortality  after  trauma. 

MATERIALS  AND  METHODS 

The  Institutional  Review  Boards  of  the  University  of 
Texas  Health  Science  Center,  Houston,  TX,  and  Brooke 
Army  Medical  Center,  Fort  Sam  Houston,  TX  approved  this 
study. 

Patients  were  identified  for  this  study  using  the  Trauma 
Vitals  database  developed  by  the  US  Army  Institute  of  Sur¬ 
gical  Research  (Fort  Sam  Houston,  TX).1  The  database  stores 
prehospital  patient  data  from  point  of  injury  until  delivery  via 
Life  Flight  helicopter  to  Memorial  Hermann  Hospital,  a  re¬ 
gional  Level  I  trauma  center  in  Houston,  TX. 

A  ProPaq  206EL  vital  signs  monitor  (WelchAllyn,  Ska- 
neateles  Falls,  NY)  and  a  standard  run  sheet  were  used  for 
data  collection.  Continuous  ECG  waveforms  were  collected 
with  an  iPAQ  (Talla-tech  RPDA,  Tallahassee,  FL)  personal 
digital  assistant  interfaced  to  the  monitor,  and  were  recorded 
at  a  sampling  frequency  of  182  Hz.  Vital  signs,  mechanism  of 
injury  (MOI,  blunt  or  penetrating),  field  Glasgow  Coma  Scale 
(GCS)  score,  Abbreviated  Injury  Scores,  Injury  Severity 
Score  (ISS),  age,  sex,  and  demographics  were  recorded. 
Blood  pressures  were  measured  automatically  by  cuff  using 
the  vital  signs  monitor.  Blood  pressures  presented  in  this 
report  were  taken  from  the  time  point  closest  to  the  ECG 
segment  extracted  for  analysis. 

Patient  Selection 

The  available  data  (117  patients  recorded  using  the 
ProPaq  206EL  vital  signs  monitor)  were  screened  for  pres¬ 
ence  of  ECG  recordings  free  of  electromechanical  noise  (se¬ 
vere  enough  to  prevent  R-wave  identification),  free  of  ectopic 
beats,  and  at  least  800  heart  beats  in  length.  Thirty-one 


patients  with  ectopy-free  ECG  waveforms  were  selected  ac¬ 
cording  to  the  above  criteria  from  the  117  candidates  avail¬ 
able.  Eighty-six  patients  were  excluded  because  of  multiple 
ectopic  beats  in  the  time  series  (S  =  34,  NonS  =  13); 
electromechanical  noise  (S  =  10,  NonS  =  22);  or  inadequate 
data  length  (S  =  5,  NonS  =  2). 

ECG  Analysis 

Eight-hundred-beat  data  sets  from  each  subject  were 
imported  into  WinCPRS  software  (Absolute  Aliens  Oy, 
Turku,  Finland)  and  analyzed  as  a  single  discrete  dataset  as 
previously  described.1518  Eight-hundred  beats  were  se¬ 
lected  because  HRV  statistics  are  affected  by  the  number  of 
data  points.  Automatic  identification  of  R  waves  was  performed 
by  the  software,  and  manually  verified  in  every  data  set.  The 
software  generated  the  instantaneous  RRI  time  series. 

The  following  are  the  main  variables  that  were  calculated: 

Time-Domain  Analyses4 

1 .  Mean  RRI:  the  mean  of  the  RRI,  measured  in  milliseconds. 

2.  RMSSD:  the  square  root  of  the  mean  of  the  sum  of  the 
squares  of  differences  between  consecutive  RRIs. 

3.  pNN50:  percentage  of  intervals  that  vary  more  then  50  ms 
from  the  previous  interval. 

Spectral  Analysis  (Fast  Fourier  Transform) 

1 .  Low-frequency  component  of  the  RRI  power  spectrum,  or 
low-frequency  (LF)  power:  influenced  by  both  sympa¬ 
thetic  and  vagal  activity. 

2.  High-frequency  (HF)  power:  influenced  by  vagal  activity. 

3.  LF/HF  and  HF/LF  ratios:  reflect  sympathovagal  balance. 

4.  Total  power:  reflects  the  strength  of  oscillations  through¬ 
out  the  entire  power  spectrum. 

5.  The  LF,  HF,  and  LF/HF  ratio  were  normalized  by  dividing 
the  LF  and  HF  spectra  by  the  total  power.  This  yielded 
normalized  powers  (LFnu,  HFnu)  and  their  ratios  (LFnu/ 
HFnu,  HFnu/LFnu). 

Spectral  Analysis  (Complex  Demodulation) 

The  method  of  complex  demodulation  (CDM)  provides 
continuous  assessment  of  the  amplitude  of  high-  and  low- 
frequency  fluctuations  in  the  RRI.19 

1.  CDM  LF:  a  measure  of  the  amplitude  of  the  low-frequency 
fluctuations  in  the  RRI. 

2.  CDM  HF:  a  measure  of  the  amplitude  of  the  high-frequency 
fluctuations  in  the  RRI. 

Complexity  Analysis  (Nonlinear  Statistics) 

1.  Approximate  entropy  (ApEn)20  and  sample  entropy 
(SampEn)21:  measure  the  amount  of  irregularity  in  the 
RRI  signal. 

2.  Fractal  dimension  by  dispersion  analysis  (FDDA)  and 
fractal  dimension  by  curve  lengths  (FDCL):  determine  the 
fractal  organization  of  the  signal.18 
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3.  Detrended  fluctuations  analysis  (DFA):  determines  fractal¬ 
like  correlation  properties  and  uncovers  short-  and  long- 
range  correlations  within  the  signal.22 

4.  Similarity  of  distributions  (SOD):  explores  the  probability 
of  similar  RR1  signal-amplitude  distributions  as  a  function 
of  time.23 

5.  Signal  stationarity  (StatAv):  assesses  whether  the  mean 
and  SD  of  the  signal  change  during  time  during  each  data 
set.24 

6.  Symbol  dynamics  indices:  symbol  dynamics  entropy 
(SymDyn),  percentage  of  forbidden  words  (FW),  and  nor¬ 
malized  symbol  dynamics  entropy  (DisnEn)  collectively 
measure  the  probability  of  certain  patterns  within  the  RRI 
time  series.25 

Statistical  Analysis 

SAS  version  8.1  (SAS  Institute,  Cary,  NC)  was  used. 
Univariate  analysis  was  performed  using  two-sample  Stu¬ 
dent’s  t  test  or  Wilcoxon’s  ranked  sum  test  as  appropriate  for 
continuous  variables,  and  the  Cochran-Mantel-Haenszel  sta¬ 
tistic  for  score  variables.  In  addition,  correlation  coefficients 
were  calculated  to  determine  relationships  between  continu¬ 
ous  variables  (Pearson  correlation)  and  between  dichotomous 
and  continuous  variables  (point-biserial  correlation).  Spear¬ 
man  correlations  were  calculated  between  ordinal  and  con¬ 
tinuous  variables. 

Multiple  logistic  regressions  with  stepwise  selection  and 
likelihood  ratio  tests  were  performed  to  identify  indepen¬ 
dent  predictors  of  mortality.  Candidate  variables  were 
ECG-derived  metrics  as  well  as  the  motor  component  of 
the  field  GCS  (GCSMOTOR)  and  the  ISS.  In  the  construc¬ 
tion  of  logistic  regression  models  to  predict  mortality,  we 
considered  this  to  represent  a  diagnostic  problem  with 
three  overlapping  phases.  In  the  first  phase  (“remote  tri¬ 
age”),  only  data  derived  from  the  RRI,  and  thus  potentially 
available  by  remote  telemetry,  were  considered.  In  the 
second  phase  (“prehospital  care”),  additional  data  avail¬ 
able  to  the  field  medic,  to  include  the  GCSMOTOR,  were 
also  considered.  In  the  third  and  final  phase  (“definitive 
care”),  data  available  during  hospitalization,  including  the 
ISS,  were  also  considered. 

We  chose  variables  with  a  p  value  of  <0.2  by  univariate 
analysis  as  candidates  for  the  logistic  models.  The  Hosmer- 
Lemeshow  goodness-of-fit  test  was  used  to  estimate  the  re¬ 
gression  model  fit.  A  receiver-operating  characteristic  curve 
was  constructed  to  assess  the  diagnostic  performance  of  pre¬ 
dictive  equations.  Estimated  odds  ratios  and  their  95%  con¬ 
fidence  intervals  (CIs)  were  determined  by  the  maximum 
likelihood  method.  The  change  in  the  Pearson  f  statistic 
caused  by  deleting  an  individual  observation  was  used  to 
detect  ill-fitted  observations  or  outliers.  If  the  model  exclud¬ 
ing  outliers  and  influential  cases  had  a  classification  accuracy 
rate  that  was  better  than  the  baseline  model,  which  included 
all  cases,  the  revised  model  was  used.  If  the  accuracy  rate  of 


the  revised  model  without  outliers  and  influential  cases  was 
less  than  2%  more  accurate,  the  baseline  model  was  retained. 

RESULTS 

Age  and  hemodynamic  data,  to  include  the  heart  rate, 
were  not  statistically  different  between  nonsurvivors 
(NonS)  and  survivors  (S)  (Table  1).  Male  gender  and  blunt 
mechanism  of  injury  were  relatively  more  common  in  the 
S  group  than  in  the  NonS  group.  NonS  were  more  severely 
injured  and  had  a  significantly  lower  GCSTOXAL  and 
GCSmotor,  and  a  higher  ISS  and  Abbreviated  Injury 
Score  for  the  head  (Table  1).  Mean  survival  time  (NonS) 
was  129  hours  ±  62  hours. 

Time-domain  and  spectral  analysis  results  are  shown  in 
Table  2.  CDM  revealed  the  amplitude  of  the  LF  oscillations 
of  the  RRI  (significantly)  and  the  HF  oscillations  of  the  RRI 
(not  significantly)  to  be  lower  in  NonS  (Table  2).  We  addi¬ 
tionally  observed  the  following  trends  (not  significant).  NonS 
had  lower  nonnormalized  HF  and  higher  HFnu  than  the  S 
group.  The  LF,  LFnu,  as  well  as  LF/HF  and  LFnu/HFnu  ratios 
were  lower  in  NonS. 

Complexity  measures  (obtained  by  nonlinear  statis¬ 
tics)  (Table  3),  to  include  ApEn,  SampEn,  FDDA  and  the 
short-term  scaling  exponent  by  DFA  were  statistically 
(p  <  0.05)  lower  in  NonS  than  in  S.  Lower  FDCL,  Lempel-Ziv 
entropy  (LZEn),  and  spectral  entropy  (SpEn)  in  NonS  were 
not  statistically  distinguishable  from  S  (data  not  shown). 
SOD  was  higher  in  NonS.  StatAv  was  not  different  between 
the  groups  (Table  3). 

Distribution  of  symbol  6  (Dis_6),  a  symbol  dynamics  mea¬ 
sure  of  RRI  complexity,  was  lower  in  NonS  (273.55  ±  27.99  vs. 
336.60  ±  16.54,  p  <  0.05).  Distribution  of  symbol  2  (Dis_2), 
was  higher  in  NonS  (292.6  ±  34.4  vs.  222  ±  21.3,  p  =  0.08). 
Other  symbol  distributions  were  not  different  between  groups 


Table  1  Demographics,  Conventional  Vital  Signs,  and 
Injury  Scores 

Variable 

NonS 

S 

Age  (yr) 

43.36  ±  5.79 

38.10  ±  3.40 

Sex 

M  (7),  F  (4) 

M  (15),  F  (5) 

MOI 

B  (6),  P  (5) 

B  (13),  P  (7) 

HR 

117.46  ±  8.54 

99.63  ±  4.39 

MAP 

74.62  ±  9.54 

82.7  ±  4.84 

GCSTotal 

8.64  ±  1.70* 

13.17  ±  0.82 

GCSMOjor 

3.36  ±  0.72t 

5.50  ±  0.32 

aishead 

3.00  ±  0.73* 

0.70  ±  0.34 

ISS 

36.55  ±  2.90* 

12.40  ±  1.97 

Data  are  means  ±  SEM. 

*p  <  0.05. 

*p  <  0.01. 

*p  <  0.001. 

NonS,  nonsurvivors;  S,  survivors;  MOI,  mechanism  of  injury;  B, 
blunt;  P,  penetrating;  MAP,  mean  arterial  pressure;  HR,  heart  rate; 
GCSTotal.  field  Glasgow  Coma  Score  total;  GCSMOTOr,  field  Glas¬ 
gow  Coma  Score  motor;  AIShead,  Abbreviated  Injury  Score  for  the 
head;  ISS,  Injury  Severity  Score. 
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Table  2  Time-  and  Frequency-Domain  Analysis  Results 


Variable 

NonS 

S 

RRI 

543.91  ±  47.37 

626.30  ±  29.47 

RMSSD 

8.09  ±  2.16 

11.60  ±  1.99 

pNN50 

0.38  ±  0.36 

1.22  ±  0.99 

TP 

1673  ± 1325 

576.35  ±  151.55 

LF 

174.64  ±  75.12 

233.00  ±  62.22 

HF 

44.45  ±  17.83 

73.05  ±  26.03 

LF/HF 

481.01  ±  201.53 

558.12  ±  175.65 

HF/LF 

0.73  ±  0.23 

0.34  ±  0.05 

LFnu 

0.57  ±  0.09 

0.74  ±  0.03 

HFnu 

0.29  ±  0.05 

0.22  ±  0.02 

LFnu/HFnu 

4.81  ±  2.01 

5.55  ±  1.73 

HFnu/LFnll 

0.67  ±  0.20 

0.33  ±  0.04 

CDM  LF 

9.09  ±  4.15* 

13.10  ±  2.11 

CDM  HF 

3.27  ±  1.05 

6.15  ±  1.24 

Data  are  means  ±  SEM. 

*p  <  0.05. 

NonS,  non-survivors;  S,  survivors; 

RRI,  mean  R-to-R  interval  of 

the  ECG;  RMSSD, 

root  mean  square 

standard  deviation;  pNN50, 

percentage  of  R-R  intervals  that  vary  by  at  least  50  ms;  TP,  total 
R-to-R  interval  spectral  power  (0.003-0.4  Hz  ms2);  LF,  RRI  spectral 
power  at  the  low  frequency  (0.05-0.15  Hz  ms2);  HF,  RRI  spectral 
power  at  the  high  frequency  (0.15-0.4  Hz  ms2);  LF/HF,  the  nonnormal- 
ized  ratio  of  LF  to  HF;  HF/LF,  the  nonnormalized  ratio  of  HF  to  LF;  LFnu 
and  HFnu,  LF  and  HF  spectral  powers  normalized  to  total  power;  LFnu/ 

HFnu,  ratio  of  LFnu 

to  HF nu;  HF nu/LF nu, 

ratio  of  HFnli  to  LFnu;  CDM 

LF,  amplitude  of  the  LF  oscillations;  CDM  HF,  amplitude  of  the  HF 

oscillations. 

Table  3  Nonlinear  Analysis  Results 

Variable 

NonS 

S 

ApEn 

0.87  ±  0.06* 

1.09  ±  0.04 

SampEn 

0.80  ±  0.08* 

1.10  ±  0.05 

FDDA 

1.08  ±  0.02f 

1.13  ±  0.01 

DFA 

0.93  ±  0.1 4f 

1.26  ±  0.08 

SOD 

0.28  ±  0.04f 

0.19  ±  0.02 

StatAv 

0.89  ±  0.05 

0.82  ±  0.03 

FW  (%) 

51.45  ±  4.91 

52.55  ±  2.28 

DisnEn 

0.62  ±  0.04 

0.65  ±  0.02 

Data  are  means  ±  SEM. 

*p  <  0.01. 

fp  <  0.05. 

ApEn,  approximate  entropy;  SampEn,  sample  entropy;  FDDA, 
fractal  dimension  by  dispersion  analysis;  DFA,  detrended  fluctuations 
analysis;  SOD,  similarity  of  distributions;  StatAv,  stationarity;  FW, 
percentage  of  forbidden  words;  DisnEn,  normalized  symbol-dynamics 
entropy. 

(data  not  shown).  FW  were  not  different  and  the  normalized 
entropy  of  the  symbol  dynamics  (DisnEn)  showed  lower  numer¬ 
ical  values  (N.S.)  for  the  NonS  (Table  3). 

Association  With  Mortality 

As  explained  above,  construction  of  logistic  regression 
models  for  prediction  of  mortality  progressed  through  three 
phases.  For  the  remote  triage  phase,  RRI  data  were  consid¬ 
ered.  ApEn  and  SampEn  were  highly  correlated  with  each 


other  ( r  =  0.99),  and  we  selected  ApEn  for  inclusion  into  our 
predictive  models.  Out  of  all  the  ECG-calculated  variables, 
two  independent  predictors  of  mortality  were  identified: 
ApEn  and  Dis_2.  The  predictive  model  based  on  these  vari¬ 
ables  produced  the  following  equation: 

P  ( mortality )  =  ek/(  1  —  ek),  where 

k  =  4.29  -  7.33  *  (ApEn)  +  0.0092  *  (Dis_2). 

The  area  under  the  receiver-operating  characteristic 
curve  (AUC)  =  0.86  (95%  Cl  =  0.71-1.0).  After  exclusion 
of  outliers  (1  case  in  NonS  and  2  in  S)  the  revised  model 
improved: 

P  ( mortality )  =  ek/(  1  —  ek),  where 

k  =  9.91  -  14.78*  (ApEn)  +  0.02  *  (Dis_2). 

AUC  =  0.956,  (95%  Cl  =  0.86-1.0;  Fig.  1). 

For  the  second  or  prehospital  care  phase,  GCSMOTOR 
was  added  to  the  model  yielding: 

P  (mortality)  =  ek/(l  —  ek),  where  k  =  5.46 

-  6.25*  (ApEn)  +  0.01  *  (Dis_2)  -  0.41  *  (GCSMOTOR). 

AUC  =  0.886  (95%  Cl  =  0.75-1.0).  After  exclusion  of 
outliers  (1  case  in  group  S)  Dis_2  was  excluded  from  the 
model,  whereas  ApEn  was  retained; 

P  (mortality)  =  ek/(l  —  ek),  where 

k  =  14.95  -  10.92  *  (ApEn)  -  0.87  *  (GCSMOTOR). 

AUC  =  0.92  (95%  Cl  =  0.80-1.0;  Fig.  2). 

Upon  consideration  of  the  ISS  for  the  third  or  definitive 
care  phase,  it  replaced  the  GCSMOTOR  in  the  model,  whereas 
ApEn  was  again  retained: 

P  (mortality)  =  ek/(  1  —  ek),  where 

k  =  -1.11  -  7.29*  (ApEn)  +  0.29  *  (ISS). 

AUC  =  0.977  (95%  Cl  =  0.91-1.0;  Fig.  3). 

For  all  the  above  models,  the  Hosmer  and  Lemeshow 
goodness-of-fit  test  revealed  no  significant  departure  from 
good  fit  (p  >  0.2). 

DISCUSSION 

This  report  introduces  the  use  of  complexity  analysis  of 
the  RRI  for  prediction  of  mortality  in  prehospital  trauma 
patients.  There  were  two  principal  findings:  (1)  prehospital 
loss  of  RRI  complexity,  as  measured  by  several  complemen¬ 
tary  but  computationally  different  nonlinear  metrics,  charac¬ 
terized  and  separated  nonsurvivors  from  survivors.  (2)  One 
measure  of  RRI  complexity,  ApEn,  was  an  independent  pre¬ 
dictor  of  in-hospital  mortality,  even  when  variables  such  as 
field  GCSmotor  and  ISS  were  taken  into  account.  ApEn 
outperformed  traditional  vital  signs,  such  as  blood  pressure  and 
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Fig.  1.  Receiver-operating  characteristic  curx’e  of  the  model  based 
on  ECG-derived  ( ApEn  and  Dis_2)  predictors  of  mortality  (“remote 
triage”  model). 
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Fig.  2.  Receiver-operating  characteristic  curve  of  the  model  based 
on  ApEn  and  GCSMOTOR  (“ prehospital  care”  model). 
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Fig.  3.  Receiver-operating  characteristic  cur\>e  of  the  model  based 
on  ApEn  and  ISS  (“definitive  care”  model). 
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heart  rate,  which  did  not  demonstrate  statistically  significant 
differences  between  nonsurvivors  and  survivors  in  this  study. 

Complexity  Analysis 

Physiologic  processes  such  as  the  RRI  display  nonlinear 
responses  to  stimuli  and  feature  complex,  irregular  patterns 


of  variability.26-28  These  complex  patterns  in  the  signal 
reflect  input  from  multiple,  interacting,  feedback-controlled 
systems.26'28'29  We  think  that  analysis  of  the  complexity  of 
the  RRI  with  nonlinear  methods  allows  assessment  of  the 
complexity,  and  thus  the  health,  of  the  underlying  system, 
which  controls  the  heart  rate. 20,27,29 

Why  does  RRI  complexity  exist?  Passage  of  information 
through  nonlinear  regulatory  networks  probably  facilitates 
adaptability  to  stress,14,28  whereas  loss  of  complexity  may 
imply  maladaptation  and  informational  isolation  of  the 
system.30  This  may  explain  why  a  system  with  higher  com¬ 
plexity  is  associated  with  survival  and  is  more  error 
tolerant,31  whereas  a  system  with  low  complexity  and  less 
variability  is  associated  with  disease13  or  death.  However,  the 
precise  mechanism(s)  responsible  for  the  observed  changes  in 
complexity  with  disease  or  shock  remain  to  be  deciphered. 

One  method  of  measuring  RRI  complexity  quantifies  the 
amount  of  irregularity,  or  entropy,  in  the  signal.  In  our  study, 
ApEn  and  SampEn  were  lower  in  nonsurvivors  than  in  sur¬ 
vivors.  Palazzolo  et  al.  documented  a  decrease  in  ApEn 
secondary  to  hypotension  in  dogs.32  Similarly,  we  showed 
decreased  ApEn  and  SampEn  during  hemorrhagic  shock  in 
anesthetized  swine.15  Hogue  et  al.  identified  that  decreased 
entropy  was  an  independent  predictor  of  atrial  fibrillation  in 
patients  after  coronary  artery  bypass  grafting.33  Decreased 
entropy  was  associated  with  myocardial  ischemia,  manifested 
by  angina.34  Thus,  decreases  in  RRI  complexity  as  measured 
by  entropy  appear  to  be  associated  with  physiologic  deterio¬ 
ration  in  a  variety  of  conditions. 

Another  method  of  measuring  RRI  complexity  takes 
advantage  of  the  fact  that  the  normal  RRI  possesses  fractal 
characteristics.  This  means  that  shorter  sections  of  the  RRI 
are  structurally  similar  to  longer  sections,  i.e.  that  the  signal 
possesses  “self-similarity”  at  both  small  and  large  scales.35 
As  measured  by  FDDA  (significantly)  and  FDCL  (nonsig- 
nificantly),  NonS  showed  a  loss  of  fractal  structure.  We  have 
previously  shown  a  decrease  in  FDDA  during  hemorrhagic 
shock  in  swine.15  West  et  al.  found  a  decline  in  RRI  fractal 
dimension  in  humans  during  central  hypovolemia  induced  by 
lower  body  negative  pressure.17  A  breakdown  in  fractal  prop¬ 
erties  may  be  an  indicator  of  a  more  regular  structural  orga¬ 
nization  of  the  signal  and  solidifies  our  findings  of  lower 
complexity  in  the  NonS  group. 

DFA  is  another  method  of  measuring  the  self-similarity 
of  fractal  processes  that  quantifies  the  short-  and  long-term 
correlations  within  the  data.22  In  this  study  we  assessed  the 
scaling  exponent  reflecting  the  short-term  correlations  in  the 
RRI  signal  and  found  it  to  be  significantly  lower  in  the  NonS, 
denoting  breakdown  of  correlations.  In  24-hour  Holter  re¬ 
cordings  of  446  survivors  of  myocardial  infarction,  Huikuri  et 
al.  found  the  DFA  short-term  exponent  to  be  an  independent 
predictor  of  death.36  Loss  of  RRI  complexity  measured  by 
decreased  short-term  fractal  scaling  properties  predicted  mor¬ 
tality  in  survivors  of  acute  myocardial  infarction.  "7  Vikman  et 
al.  found  that  a  decrease  in  RRI  ApEn  and  short-term  RRI 
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dynamics  (DFA)  preceded  the  spontaneous  onset  of  atrial 
fibrillation  in  patients  with  no  structural  heart  disease, 
whereas  traditional  HRV  metrics  did  not  change.38  Our  find¬ 
ings  are  therefore  consistent  with  the  literature  indicating  an 
association  between  loss  of  RRI  fractal  scaling  and  life- 
threatening  disease. 

SOD  is  a  method  that  allows  for  rapid  online  analysis 
of  small  datasets,  as  it  explores  the  probability  of  similar 
RRI  signal  amplitude  distributions  as  a  function  of  time.23 
SOD  was  higher  in  NonS,  reflecting  greater  regularity  of 
signal  distribution  and  thus  a  state  of  lower  complexity 
among  NonS. 

The  symbol-dynamic  indices  collectively  measure  the 
probability  of  certain  patterns  within  the  RRI  time  series.25  In 
our  study,  the  FW  and  the  normalized  entropy  of  the  symbol 
dynamics  (DisnEn)  were  not  different  between  the  groups, 
possibly  denoting  low  sensitivity  of  this  method.32  In  addi¬ 
tion,  we  found  that  the  average  of  pattern  6  of  symbol  words 
(Dis_6)  was  higher  and  pattern  2  (Dis_2)  was  lower  in  NonS, 
with  the  latter  included  into  the  “remote  triage”  predictive 
model.  At  this  point,  limited  experience  with  these  methods 
prevents  us  from  drawing  firm  conclusions  about  the  physi¬ 
ologic  meaning  of  these  specific  patterns  (Dis_6  and  Dis_2), 
and  application  of  the  methods  to  larger  data  sets  will  be 
needed  to  place  these  findings  in  perspective. 

Spectral  Analysis 

The  method  of  CDM  provides  continuous  assessment  of 
the  amplitude  of  high-  and  low-frequency  fluctuations  in  the 
RRI,  and  is  insensitive  to  data  length  and  nonstationarity  in 
the  signal.19  As  measured  by  CDM,  the  amplitude  of  LF  and 
HF  oscillations  (the  latter  nonsignificantly)  were  both  lower 
in  NonS.  These  changes  paralleled  those  seen  in  the  standard 
nonnormalized  frequency-domain  metrics.  Overall,  these 
findings  indicate  inappropriately  low  sympathetic  activity  in 
the  NonS  group. 

Previous  work  by  Winchell  and  Hoyt  on  intensive  care 
unit  patients,9  as  well  as  from  our  institution  on  prehospital 
trauma  patients,10’11  identified  lack  of  sympathetic  tone  and  a 
state  of  decreased  sympathovagal  balance  in  NonS.9’10  Con¬ 
sistent  with  those  reports,  in  this  study  the  trends  in  the  LF, 
LFnu,  and  LF/HF  ratio  were  all  numerically  lower  in  the 
NonS  group  along  with  a  higher  HFnu,  HF/LF  ratio,  and 
HFnu/LFnu  ratio  (Table  2).  Lack  of  statistical  significance  in 
these  observations  likely  reflects  sample  size,  dataset  length 
and,  perhaps,  the  methodologic  limitations  of  these  methods. 

Prediction  of  Mortality 

The  first  predictive  model  explored  in  our  study  reflected 
the  discriminative  capacity  of  ECG-derived  metrics  alone, 
and  modeled  a  scenario  in  which  an  ECG  signal  is  analyzed 
remotely  and  before  hands-on  evaluation  of  the  casualty  (“re¬ 
mote  triage”  model).  The  respectable  sensitivity  and  speci¬ 
ficity  (that  further  improved  with  elimination  of  outliers) 
underlines  the  potential  utility  of  entropy  measurements  for 


remote  prognosis  and  diagnosis.  Inclusion  of  the  GCSMOTOR 
in  the  second  predictive  model  conceptualizes  arrival  of  a 
medic  to  the  casualty  and  a  physical  examination  in  addition 
to  the  ECG  analysis  (“prehospital  care”  model).  Finally,  in 
the  third  (“definitive  care”)  model,  ApEn  was  retained  and 
ISS  was  added  as  independent  predictors  of  mortality.  Thus, 
in  this  dataset  ApEn  possibly  accounts  for  differences  in  the 
robustness  of  an  individual’s  response  to  injury,  even  when 
ISS  is  taken  into  account. 

Limitations 

The  hemodynamic  data  (heart  rate  and  mean  arterial 
blood  pressure)  were  numerically  but  not  statistically  differ¬ 
ent  between  S  and  NonS  in  this  study.  Thus,  the  possibility  of 
a  type  II  error  exists  in  this  small  sample  size.  We  investi¬ 
gated  prehospital  trauma  patients  with  ECGs  that  did  not 
contain  prohibitive  levels  of  mechanical  noise  or  arrhythmias, 
and  were  800  beats  long.  These  limitations  led  us  to  the 
sample  size  of  3 1  of  the  original  1 1 7  as  outlined  in  our  patient 
inclusion  criteria.  It  is  possible  that  those  trauma  victims  who 
are  more  critically  ill,  and  therefore  more  likely  to  die,  might 
have  more  noisy  ECGs  because  of,  for  example,  ongoing 
interventions  during  transport.  Indeed,  the  mortality  rate  was 
43%  in  the  excluded  group  versus  35%  in  the  included  group. 
Otherwise,  excluded  patients  did  not  differ  from  included 
patients  with  respect  to  the  male-to-female  ratio,  age,  mean 
arterial  pressure,  or  ISS.  Clearly,  however,  it  would  be  pre¬ 
mature  to  generalize  our  findings  to  other  trauma  patients. 
Efforts  are  ongoing  to  improve  the  quality  of  ECG  results 
collected  in  the  field,  and  to  determine  the  absolute  minimum 
dataset  length  (number  of  beats)  needed  for  accurate  patient 
classification.  Reduction  of  dataset  size  would  widen  the 
utility  of  the  proposed  methods  during  short  transport  times. 

The  endpoint,  in-hospital  mortality,  was  chosen  in  our 
study  for  its  unequivocal  ability  to  categorize  outcome.  How¬ 
ever,  of  greater  importance  to  the  prehospital  provider  is  the 
ability  to  predict  the  need  for  a  lifesaving  intervention.  Stud¬ 
ies  are  ongoing  to  verify  the  utility  of  the  methods  explored 
here  in  prediction  of  lifesaving  interventions. 

The  analysis  methods  used  in  this  study  require  modest 
computing  power  and  could  be  incorporated  into  microchips, 
personal  digital  assistants,  or  commercial  monitors,  making 
these  tools  readily  available  for  clinical  use.  It  is  clear  that 
there  are  several  quite  different  ways  to  quantify  what  is 
generally  called  “heart-rate  variability”,  to  include  indices 
derived  from  spectral  analysis  (such  as  sympathovagal  bal¬ 
ance)  and  from  complexity  analysis  (such  as  approximate 
entropy).  Prospective,  large-scale  clinical  trials  will  be 
needed  to  determine  the  clinical  utility  of  these  methods  in 
guiding  decision  making,  and  to  define  the  merits  of  these 
different  approaches. 

CONCLUSION 

Loss  of  RRI  complexity,  as  evidenced  by  decreased 
entropy,  diminished  fractal  scaling,  loss  of  short-term  cor- 
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relations  in  the  RRI,  and  higher  similarity  of  distributions 
in  the  signal,  characterized  nonsurviving  trauma  patients 
and  distinguished  them  from  survivors.  The  fact  that  sev¬ 
eral  computationally  different  methods  pointed  to  loss  of 
complexity  in  NonS  lends  additional  weight  to  these  find¬ 
ings.  Prehospital  loss  of  RRI  complexity,  measured  by 
ApEn,  was  associated  with  mortality  independent  of  the 
field  GCSmotor  and  hospital  ISS.  Evaluation  of  cardio¬ 
vascular  complexity  may  be  useful  for  noninvasive,  remote 
triage. 
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